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Abstract 

Because of the quantum fluid properties of a neutron star core's neutrons and pro- 
tons, its magnetic field is expected to be coupled strongly to its spin. This predicts 
a simple evolution of the surface- field of such stars as they spin down or, less com- 
monly, are spun up. Consequences and comparisons with observations are given for 
properties of solitary spinning down pulsars, including their glitches and spin-down 
ages, X-ray pulsars, and the formation and pulse characteristics of Millisecond Pul- 
sars. For none of these is there a present conflict between model predictions and 
what has been observed. 
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1 Introduction 



Enough is understood about the dynamics of the components of a standard 
(non- magnet ar, non-strange) neutron star (NS) to support what should be a 
reliable description of what happens within a spinning magnetized NS as it 
ages and spins down or, in rarer cases, when it is spun up. 



In a cool core below the crust of a spinning NS superconducting protons coexist 
with more abundant superfluid neutrons (SF-n) to form a giant atomic nucleus 
which contains within it a neutralizing sea of relativistic degenerate electrons. 
The neutrons rotate with a spin-period P (sec) = 2tc/Q only by forming 
a nearly uniform array of corotating quantized vortex lines parallel to the 
spin axis, with an area density n v ~ 10 4 cm -2 P _1 . The array must contract 
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Fig. 1. A moving quantized vortex- line in a NS core's superfluid neutrons pulling 
a pair of the core's proton superfluid quantized flux-tubes anchored in the star's 
solid, conducting crust (shown dotted). 

(expand) when the NS spins up (down). In stellar core neutron spin-up or spin- 
down, a vortex a distance r± from the spin-axis generally moves outward with 
a velocity v v = r±(P/2P) until r± reaches the core's neutron superfluid radius 
(R). Any stellar magnetic field passing below the stellar crust must, in order to 
penetrate through the core's superconducting protons (SC-p), become a very 
dense array of quantized flux-tubes (n$ ~ 5 x 10 18 B\2 cmT 2 with B the local 
average magnetic field). Each tube carries a flux 2 x 1CT 7 G cm 2 and a magnetic 
field B c ~ 10 15 G. 1 The initial magnetic field within the core of a neutron star 
is expected to have both toroidal and very non-uniform poloidal components. 
The web of flux-tubes formed after the transition to superconductivity is then 
much more complicated and irregular than the neutron vortex-array as well as 
of order 10 14 times more dense. Because of the velocity dependence of the short 
range nuclear force between neutrons and protons, there is a strong interaction 
between the neutron-superfluid's vortex-lines and the proton-superconductor's 
flux-tubes if they come closer to each other than about 10~ 11 cm. Consequently, 
when Pj£ flux tubes will be pus hed (or pulled) by the moving neutron 
vortices@; HBSBSQBEO)- A realistic flux-tube array will be 
forced to move along with a changing SF-n vortex array which threads it as 
long as the force at a vortex-line flux-tube juncture does not grow so large that 



This assumes Type II proton superconductivity in the NS core, below the crust, 
the common result of many calculations. If it were Type I, with many thin regions 
of B > several B c , and B ~ in between (flih. the impact on surface B of changing 
NS spin proposed below would not change significantly. If, however, (B), the locally 
averaged B inside the NS core exceeds a critical field somewhat greater than B c , 
the core's protons would not become superconducting. This may well be the case 
for most (or all) "magnetars" (jj^ ) . 



2 



vortex lines cut through flux-tubes. The drag on moving flux-tube arrays from 
their small average velocities (f± < 10 _5 cm s _1 ) in spinning-down pulsars, 
cool (old) enough to have SF-n cores, seems far too small to cause such cut- 
through. 2 

The main quantitative uncertainty in the model described below is the maxi- 
mum sustainable shear-strain (8 m ~ 10~ 4 — 1CT 3 ?) on the conducting crust, 
which anchors core flux-tubes (cf. Fig 1), before the crust yield-strength is 
exceeded. An estimate (Jiih for that maximum sustainable crustal shear-stress, 
compared to that from the (B 2 )/8tt ~ (B)B c /8tc of the core's flux-tube array, 
supports a NS model in which the crust yields before the core's flux-tubes are 
cut through by its moving SF-n vortex array, as long as Byi > 1. Even for 
much smaller £>i 2 , flux-tube anchoring by the conducting crust would result 
in such cut-through only when the NS's spin-down age (P/2P) exceeds the 
crust's Eddy current dissipation time (~ 10 7 yrs.). Then in most observation- 
ally relevant regimes the motion of the magnetic flux-tube array near the top 
of the NS core (and B at the crust surface above it) follows that of the SF-n 
vortex array which threads it. This forms the basis of a very simple model 
for describing predicted changes in pulsar magnetic fields during NS spin-up 
or spin-down which agrees well with a variety of different families of pulsar 
observations. 



2 Magnetic field changes in spinning up neutron stars 



NS spin-up, when sustained long enough so that one of the above criteria 
for relaxation of shear-stress from crust-anchored magnetic flux before cut- 
through is met, leads to a "squeezing" of surface B toward the NS spin- axis. 
After a large decrease in spin-period from an initial Pq to P <C Pq all flux 
would enter and leave the core's surface from the small area within a radius 
R(P / Pq) 1 / 2 of the NS's spin-axis. This surface B-field change is represented in 
Figs 2-3 for the special case when the magnetic flux which exits the NS surface 
from its upper (lower) spin- hemisphere returns to the stellar surface in it's 
lower (upper) one. Potentially observable features of such a "spin-squeezed" 
surface B configuration include the following. 

a) A dipole moment nearly aligned along the NS spin-axis. 



2 Jones (UJ) has recently found that electron scattering on flux-tube cores allows 
easier passage of flux-tubes through the SC-p than had been previously estimated 
(e.g. (25)). In addition, an expected motion- induced flux-tube bunching instability 
would allow easy co-motion of flux-tubes with the local electron plus SC-p fluid in 
which they are embedded |26h. 
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Fig. 2. A single flux-tube (one of 10 31 ) Fig. 3. The flux-tube and vortex array 
and some of the NS's arrayed vortices of Fig. 2 after a large stellar spin-up. 
(8 of 10 17 ). 



b) A greatly diminished polar cap spanning the "open" field lines when 
P/P -> 0. For P < P (A/R)(tt R/c) 1/2 with A the crust thickness 
(~ lO^R), the canonical polar cap radius, r p = R{VtR/ c) 1 ^ 2 , shrinks 
to r' p = A(QR/cy/ 2 

c) A 5-field just above the polar cap which has almost no curvature. 

If the pre-spin-up surface B has a sunspot-like configuration (i.e. flux returning 
to the NS surface in the same hemisphere as that from which it left), the spin- 
up-squeezed field change is represented in Figs 4 and 5. In this case, potentially 
observable features when P <C P include the following. 

d) A pulsar dipole moment nearly orthogonal to the NS spin-axis, and 

e) positioned at the crust-core interface. 

f) A dipole moment (//), or more precisely the component of fi per- 
pendicular to ft, reduced from its pre-spin-up size: 

fi±(Po) ~ \Po) 

A more general (and very probably more realistic) pre-spin-up configuration 
has flux emitted from one spin-hemisphere returning to the stellar surface in 
both, as in Fig. 6. Spin-up squeezing then typically gives the surface field 
configuration represented in Fig. 7, a spin-squeezed, nearly orthogonal dipole 



(1) 



4 




Fig. 4. A single flux-tube, part of a 
sunspot-like -B-field geometry in which 
flux from a spin-hemisphere of the sur- 
face returns to the surface in that same 
hemisphere. 




Fig. 5. The flux-tube and vortex array 
of Fig. 4 after a large stellar spin-up. 



on the NS spin-axis with properties (d), (e), and (f), together with an aligned 
dipole on the spin-axis whose external field is well-represented by North and 
South poles a distance 2(R — A) apart. Further spin-up could lead to the Figs. 
8 and 3 configuration; that of Fig. 9 and 5 would be realized only if S2 of Figs. 
6 and 7 is negligible. 
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Fig. 8. The field from Fig. 7 after further Fig. 9. The field from Fig. 6 after large 
spin-up. spin-up when the S2 contribution to 

Fig. 7 is negligible. 



3 Magnetic field changes in spinning down neutron stars 

Consequences of the coupling between a spin-down expansion of a NS's SF- 
n vortex-array and its SC-p flux-tubes should appear in several observable 
phases which begin after the NS has cooled enough that the vortex-line array 
and the flux-tube one have both been formed (typically after about 10 3 yrs.). 

a) As in Eqn (1), except that P > P , ^±(P) initially grows as P 1 / 2 . 
This increase is initially not sensitive to the configuration of surface 



b) When P ~ several Po, a good fraction of a NS's core flux-tubes will 
have been pushed outwards from the spin-axis to r± ~ R. These 
cannot, of course, continue to move outward (Fig. 11) so that Eqn 
(1) no longer holds. Rather, the mixture of expanding and crust- 
constrained flux-tubes gives: 



with the exact value of n dependent on details of a core's 5-field 
configuration. 

c) The crust can delay, but not indefinitely prevent, expulsion of this 



B (cf. Fig. 10). 




(0 < n < 1/2) 



(2) 
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Fig. 10. The flux-tube and vortex-array 
of Fig. 4 after some spin-down. The ex- 
panded configuration would not differ in 
an important way if it had begun as that 
of Fig. 2. 



Fig. 11. The equatorial plane, viewed 
from above, of a configuration like that 
of Fig. 10, but with two flux-tubes. One 
tube is being expelled into the crust by 
the expanded vortex array and will ul- 
timately be eliminated by reconnection. 



magnetic field from the NS. When P ~ several P , intertwined vor- 
tex plus flux which have been pushed into the core-crust interface 
will stress the crust enough to exceed its shear-strength (Sect. 4 and 
Figs. 15 and 16). Then crust movements begin that lead to i?-field 
reconnections. Flux that is threaded by SF-n vortex lines that have 
not yet reached rj_ ~ R, and thus have not yet disappeared in this 
way, are the remaining source for the NS's dipole moment. The sum 
of all this remaining flux oc the total number of remaining vortex- 
lines (oc Q). Then, Eqn (2) holds with n — —l. 

d) When this remaining B at the crust bottom (oc Q) drops to and be- 
low ~ 10 12 G, shear-stress on the crust would no longer be expected 
to exceed the crust's yield-strength. The NS's surface B may then 
lag that at the base of its crust by as much as 10 7 yrs., the crust's 
Eddy current dissipation time. 



4 Comparisons of pulsar expectations with model expectations 

Fig. 12 shows observationally inferred surface dipole fields (B) as a function 
of their P for about 10 3 radiopulsars (B is calculated from measured P and 
P, IQ = -ni VI 3 c" 3 ; B = fi ± R~ 3 and I = 10 45 g cm 2 .). Segments of B(P) 
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Fig. 12. Dipole-1? observed on pulsar surfaces (inferred from measured P and P) as 
a function of pulsar period (P) (|12l). The solid line segments are the evolutionary 
segments for B of Sect. 4, based upon the model of Sects. 2 and 3. The dash-dot 
diagonal is the steady state spin-up line from an accretion disk. The horizontal 
(D — > D') is that for a NS surface above a core surface {D — > i?). 

based upon the model of Sects 2 and 3, are shown for a typical pulsar by the 
doubled or single solid lines. 

(1) Point A is the (B,P) where, typically, flux-tubes and vortex lines begin 
coexistence. 



(2) (A — > C) is the expanding array prediction of Sect. 3(b): B oc P n with the 
model prediction < n < 0.5. The index n is known o nly in the several 
cases where P is also measured: n = +0.3,0.1,0.6,0.1 QUI Hi). 



(3) (C — > D) is the flux-expulsion and continual reconnection segment of 
Sect. 3(e). The model predicts (n) = — 1 for h averaged over the (C — > D) 
history of any one pulsar. Reliable P are not generally measurable in this 
(B,P) region. However comparison of the spin-down times P/2P with 
actual pulsar ages, inferred from probable distance traveled since birthjl), 
give (n) = —0.8 ± 0.4, not inconsistent with the model prediction. 



(4) (D — > E) is the core-surface/crust-base B evolution for ~ 10 10 yrs. The 
horizontal (D — > D') is the NS crust's surface field, remaining near 10 12 
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G for ~ 10 7 yrs. as discussed in Sect. 3(d). This segment should be 
characteristic of typical "X-ray pulsars" (NSs in binaries spun up or down 
by active companions through a wide range of P (e.g. Hercules X-l with 
P ~ Is to Vela X-l with P ~ 10 3 s) until crustal Eddy current decay 
allows a (D' —>■ E) decay from some D' region. 

A small minority of NSs, after (D' — > E) segments, will be resurrected 
by accretion from a previously passive White Dwarf companion which 
now overflows its Roche lobe (LMXBs). These NSs have entered into the 
spin-up phase of Sect. 2 until they reach a steady state on the canonical 
"spin-up line" represented by the dot-dashed diagonal of Fig. 12 (for 

M = 10 _1 MEddington) ■ 



(5) {E — > F — > H) is the spin-up segment when the NS surface B has the 
sunspot geometry of Figs. 4, 5, and 9, which allows spin-up to minimal P 
before spin-up equilibrium is reached. Observations of maximally spun- 
up millisecond pulsars (MSPs) support the Sect. 2 model for such MSP 
formation: Sect. 2(d) 's high fraction of MSPs with two subpulses 180° 
apart, characteristic of orthogonal rotatorsp; 0; 11 ; BO); Sect. 2(e) 's 



B-field geometry, from linear polarization and its frequency dependence 
in such subpulses (j~ 



(6) (E — * F — > K) is the track of surface B (here the total dipole field) pre- 
dicted after large spin-up from (E) with Fig. 6 geometry to (F) with Fig. 
7 geometry. Further spin-up diminishes only the orthogonal component 
of /J, until an almost aligned rotator (Figs. 3 and 8) results when (K) is 
reached. X-ray emission from the almost aligned MSP PSR 0437 (P = 6 
ms) supports a (predicted) tiny polar cap area about (A/R) 2 ~ 1CT 2 that 
from a central dipole moment configuration for the same P (Sect. 2(b) 
and refs Q; Q)). 

Expected consequences for pulsar diplole-B changing according to the 
Sects. 2-3 model and Fig. 12 are supported by many kinds of observations. 
However, for almost all there is usually another popular explanation (e.g. 
B getting from (D) to (H) just by burial of B by accreted matter from 
a companion (ji^ Efi)). 



5 Pulsar spin-period glitches from spin-induced Z?-field changes 



Moving core flux-tubes continually build up shearing stress in the conducting 
crust which anchors 5-field that traverses it. If this stess grows to exceed the 
crust's yield strength, subsequent relaxation may, at least partly, be through 
relatively sudden crustal readjustments ("crust-breaking"). Such events would 



9 



GLITCHES 



10' 



10" 



a 1 °" 



10" 



JO 




10 1D" 
Pulsar Period (s) 



10" 



Fig. 13. Observed jumps ("glitches") in pulsar spin-rates (Afi/fl) of pulsars with 
various periods (P). The Vela- like family (V) has Afi/fl ~ 10~ 6 . The Crab- like one 
(C) has Afi/O ~ 1(T 7 - 10" 8 (0;[lB;S;0;|3). 



cause very small spin-up jumps in spinning-down NSs (spin-period "glitches"). 
The Sect. 2-3 model for the evolution of a core's flux-tube array suggests glitch 
details in pulsars similar to those of the two observed glitch families: Crab-like 
glitches (C) and the very much larger giant Vela-like ones (V) of Fig. 13. 

a) Crab-like glitches In both the (A — > C) and (C — > D) segments of 
Fig. 12, an expanding quasi-uniform vortex-array carries a flux-tube 
array outward with it. If growing flux-tube-induced stress on the 
crust is partly relaxed by "sudden" outward crust movements (of 
magnitude s) where the stress is strongest (with density preserving 
backflow elsewhere in the stratified crust) the following consequences 
are expected: 



(1) a "sudden" permanent increase in fi±, spin-down torque, and 
|f2| : Att/tt ~ s/R ~ AO (strain relaxation) < 6> max ~ 10~ 3 . 
(This is the largest non-transient fractional change in any of 
the pulsar observables expected from "breaking" the crust.) A 
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Fig. 14. The difference between Crab pulsar periods observed over a 23 yr interval 
and those predicted from extrapolation from measurement of P, P, and P at the 
beginning of that interval. These "sudden" permanent fractional jumps in spin-down 
rate (Att/n ~ +5 x 1(T 4 ) occur at glitches (AQ/Q ~ 1(T 8 - 1CT 7 ) but are 10 4 



times greater in magnitude(|la: ll 7T) . 



permanent glitch-associated jump in NS spin-down rate of this 
sign and magnitude (~ 3 x lO -4 ) is indeed observed in the larger 
Crab glitches (Fig. 14) B H B 0). 



(2) a "sudden" reduction in shear stress on the crust by the flux- 
tubes attached to it from below. This is matched by an equiv- 
alent reduction in pull-back on the core's expanding vortex ar- 
ray by the core flux-tube array attached to it. The n- vortices 
therefore "suddenly" move out to a new equilibrium position 
where the Magnus force on them is reduced by just this amount. 
The high density SF-n sea therefore spins down a bit. All the 
(less dense) charged componentes of the NS (crust, core-p and- 
e) together with the flux-attached n-vortex-array spin-up much 
more. (The total angular momentum of the NS does not change 
significantly in the brief time for development of the glitch.) A 
new equilibrium is established in which the charged components 
(all that can be observed, of course, is P of the crust's sur- 
face) have been spun up. For Crab B and P, the estimated (jl^) 
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AQ/Q ~ 10~ A (AQ/Q), consistent with both the relatively large 
Crab glitches of Fig. 14 and also with much smaller Crab glitches 
not shown there (34). 




Fig. 15. The configuration (top view) Fig _ 16 A side view of the re p resent a- 

of Fig. 11 after further spin-down. tion of the Fig> 15 confi g ura tion with 

Flux-tubes are piling up in an equatorial the addition of one furx-tube, which 

annulus at the core-crust interface. The ^he expanding vortex- array has not yet 

blocked flux-tubes, in turn, block short forced out to a radius ~ R 
segments of vortex lines which forced 
them into this annulus. 



b) Giant Vela-like (V) glitches. The second (V)-family of glitches differs 
from that of Crab-like ones (C) in several ways. 



;i) (aq/q) v ~ io 2 x (Afl/n) 



c- 



(2) V-glitches develop their AQ in less than 10 2 sec: the Afl of a V- 
glitch is already decreasing in magnitude when first resolved (jlfih . 
while C-glitches are still rising toward their full AQ for almost 



10 5 sec()7|;|2lj). 



(3) V-glitches are observed in pulsars (mainly, but not always) in 
Fig. 12 along (C — > D) while C-glitches are observed all along 
(A^C -> D). 

(4) The C-glitch proportionality between AQ/Q and AQ/Q would 
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greatly overestimate (Af2/f2) for V-glitches. 




Fig. 17. A schematic representation of 
a young NS's magnetic field just be- 
fore the NS cools to the transition tem- 
perature for proton superconductivity. 
Some shearing stress preventing an even 
more stabilized configuration is proba- 
bly borne by the NS crust which solidi- 
fied much earlier. 



Fig. 18. A representation of the Fig. 
17 magnetic field after core flux-tube 
formation and relaxation to a new 
quasi-equilibrium. The initially in- 
creased stress in the crust (cf. Fig. 
1) is assumed to exceed the crust's 
shear-stress yield strength. A later for- 
mation of SF-n vortex-lines would halt 
such relaxation. 



The existence of a second glitch family, with V-properties, is 
expected from a second effect of vortex-driven flux-tube move- 
ment in a NS core. If there were no very dense, comoving, flux- 
tube environment around them, outward moving core-vortices 
could smoothly shorten and then disappear entirely as they 
reached the core's surface at its spin-equator. (We ignore crustal 
SF-n here.) However, the strongly conducting crust there resists 
entry of the flux-tubes which the vortices also bring with them to 
the crust's base. This causes a pile-up of pushed flux-tubes into 
a small equatorial annulus (Figs. 15 and 16) which delays the 
final vortex-line disappearance. The vortex movement in which 
they vanish occurs either in vortex-line flux-tube cut-through 
events, or, more likely, in a sudden breaking of the crust which 
has been overstressed by the increasing shear-stress on it from 
the growi ng a nnulus. Giant V-glitches were proposed as such 
events ljifjl : liih . allowing a "sudden" reduction of part of this 
otherwise growing annulus of excess angular momentum and 
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also some of the magnetic flux trapped within it. These would 
not begin until enough vortex-lines, initially distributed almost 
uniformly throughout the core, have piled up in the annulus 
for the flux-tubes they bring with them to supply the needed 
shear stress. Estimates of V-glitch AQ/Q magnitudes are less 
reliable than those for C-glitch ones. A very rough one, based 
upon plausible guesses and an assumed Q/R about the same as 
those in the larger C-glitches, suggest V-glitch repetition rates 
and magnitudes not unsimilar to observed onesljia: Eil). 



ioun'flMe'fif'Xl 



Iff 



10 







to* 



*}E 1701 



t 



ID* 10* 
t*«* (***) 



Fig. 19. Observed spin-down times for pulsars (P/2P) vs the time since birth of 
these same pulsars as inferred by the ages of the supernova remnants in which they 
are still embedded fejvij)@;S- 
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6 In the beginning 

The proposed spin-down biography of a NS surface B presented in Sects. 
3,4, and 5 began at (A) (or perhaps A') in Fig. 12 when that typical NS is 
expected to be about 10 3 yrs old. Before that its crust had solidified (age 
~ a minute), its core protons had become superconducting (~ 1 yr?), and 
core neutrons became superfluid (~ 10 3 yrs?). If so, there would be a nearly 
10 3 year interval between formation of the NS core's magnetic flux-tube array 
and control of that array's movement by that of a SF-n vortex array. During 
that interval an early magneto-hydrodynamic equilibrium involving poloidal 
and toroidal fields, and some crustal shear stress (Fig. 17) would be upset 
by the dramatically altered 5-field stresses after flux-tube formation (1281) . The 
subsequent jump in shearing stress on the crust surface B change. The recent 
reconsideration of drag on moving flux-tubes (fl2l) suggests the core flux-tube 
adjustment can take ~ 10 3 yrs. For many NSs, depending on historical details 
of their B structure, dipole moments should become much smaller (Fig. 18). 
Their post-partem values and subsequent expected drops in their sizes have 
been estimated and proposed (28) as the reason many young pulsars have 
spin-down ages (P/2P) up to 10 2 times greater than their true ages (Fig. 19). 
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